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The problem of the observable equilibrium domain structure in pure antiferromagnets (and 
other thermoelastics) is investigated with the use of continuous elasticity theory. It is shown that 
completely rigid surface produces the imaginary "incompatibility elastic charges" analogous to the 
surface magnetic charges in ferromagnets. Corresponding long-range field is shown to contribute 
into the "stray" energy of the sample that governs an appearance of the domain structure, the 
contribution from the "elastic charges" being proportional to the sample volume. Competition 
between the elastic "stray" field that favors inhomogeneous strain distribution, and external field 
that tends to homogenize the sample, provides the reversible reconstruction of the domain structure 
under the action of external magnetic field. 
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I. INTRODUCTION 

The problem of antiferromagnetic (AFM) domains has 
a very long history and was discussed for many times 
(see, e.g., Rcf. |l]-fl|). This problem is tightly related with 
the question of the domain structure (DS) origin in the 
substances that suffer thermoelastic phase transitions, 
i.e. transitions at which the principal order parameter 
is symmetrically conjugated with the shear components 
of the strain tensor. The transitions of such a kind are 
usually accompanied by appearance of the spontaneous 
strains along with the simultaneous change of the elec- 
tronic structure (martensitic phase transitions) or onset 
of the vector or tensor order parameter of other nature 
(like magnetization, polarization, etc.). In the case when 
the values of the principal order parameter and/or corre- 
sponding strain tensor of low-symmetry (i.e., magnetic, 
ferrielectric, martensitic) phase are degenerated, the DS 
is thermodynamically equilibrium and reversible under 
the action of the external fields that can result in homoge- 
nization of the sample. The nature of the equilibrium DS 
of ferro- and ferrimagnets with opposite direction of mag- 
netization inside the neighboring domains is well studied 
and is usually related with the finite size of the sam- 
ple and demagnetization effect resulting from the surface 
magnetic charges. So, thegCtuestion is "What is the rea- 
son of the clearly observedETEl equilibrium DS in the cases 
when the noncompensated magnetization is equal to zero 
or vanishingly small (e.g. in weak ferromagnets)?" In 
the Ref. || we have supposed that one of the possible 
mechanism of equilibrium DS formation can be due to 
the surface effects, namely, the surface energy decrease 
yielded by the global variation of the sample shape. Al- 
though such a model has given qualitative and quantita- 
tive agreement with the magnetostriction experimentsEl, 
it application, as it is evident from theoretical point of 
view, is restricted to the rather small sanaales for which 
the surface effects can be clearly detectecO. 

In the present paner we make an attempt to generalize 
our previous resultsu and to show that the formation of 
equilibrium DS in the course of thermoelastic phase tran- 
sition can be explained by the effect of low-deformable 
sample surface, and what is principal, corresponding con- 
tribution into the sample energy being proportional not 
to the surface area but to the sample volume. 



II. MODEL 

Let us consider a sample with the volume V confined 
with the surface S. Suppose, that above the critical tem- 
perature the sample is non-stressed and non-strained. 
Below the critical temperature the crystal lattice be- 
comes unstable with respect to appearance of the spon- 
taneous strains related with the fluctuations of the order 
parameter (e.g., magnetizatioon, polarization, etc). The 
spontaneous strains can result from the nonlinear elastic 



effects, as in the case of martensites, as well as be related 
with the appearance of spontaneous AFM moment, both 
cases allow no demagnetization effects. 

In the infinite sample the only equation that defines 
distribution of spontaneous strain tensor u(r) (r is a 
space coordinate) arises from the condition of the elastic 
free energy minimum with respect to u: 



dF el (u) 
du(r) 



= 0, 



(1) 



where F e \(u) is the free energy density. Equation (jl]) 
is equivalent to the condition of the absence of internal 
forces: divcr(r) = 0, where a is a stress tensor. Elastic 
free energy density F e \(u) is usually independent upon 
the space derivatives of the strain tensor; so, the equi- 
librium strain is homogeneous, and thus trivially satisfies 
the compatibility conditionsLH that in approximation of 
linear elasticity takes a form: 



inc u(r) = — rotrot u(r) = 0. 



(2) 



If we take into account the surface of the sample, given 
parametrically by equation r = rg, the situation changes 
in a crucial way. The elastic properties of the sample 
surface differ from that of the bulk, and homogeneous 
strain defined from inevitably brings about the addi- 
tional stresses in the sample. In formal, the surface can 
be accounted for through the surface energy ternJlj: 



dS 



a s +j3 s u(r s ) 



(3) 



where as is a surface tension coefficient, (3s is the tensor 
of the surface elastic modula, and u(ts) is a surface strain 
tensor. It should be emphasized that the first term in (j^) 
is responsible for the set up of the equilibrium shape of 
the sample (it establishes the minimum surface area at 
given volume), the second term can be important when 
considering thermoelastic phase transition. 

In what follows we will suppose that the shape of the 
sample is given technologically and does not relax during 
the time of experiment, and so, will disregard an influence 
of the surface tension (||). In other words, we consider the 
energy barrier for shape variation in the strained sample 
as high enough, so, that corresponding relaxation time is 
much greater than the inverse energies characteristic to 
the phase transition. This is the case for AFM transi- 
tions at low temperature and also for martensitic phase 
transitions that usually proceed with the rate close to 
sound velocity. 

In the vicinity of the surface the spontaneous strains 
that occur during the thermoelastic phase transition can 
be represented as a sum of three summands: 

U(r) = ^ITrU + U sh car 

+ {(n s u) <g> n s + n s <g) («n s ) - n s <g> n s (n s un s )} , (4) 
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where 1 is a unit matrix, rig is a surface normal. The first 
term in (^) is an isotropic volume striction, the last three 
terms in brackets form the component that produces the 
shift along the surface normal ns, and the second, trace- 
less, term corresponds to the shear strains (with respect 
to rig). The first component is significant at the mag- 
netic phase transitions and can be neglected during the 
martensitic ones, but at any rate it does not contribute 
into the symmetry change of the crystal and hence, as 
well as the last term, could not be removed by onset of 
DS in the sample. So, in what follows, we consider only 
the shear part of the strain tensorE3 u s h C ar and define the 
strain state with respect to isomorphous magnetostric- 
tion which contribution can be taken into account triv- 
ially. 

If one assumes that the sample surface is not deformed 
during the phase transition, i.e., (3s = 0, we should im- 
pose the standard boundary conditions for zero external 
stress 



\ du 



n s (rs) = 0, 



(5) 



res 



and additional conditions of the absence of surface 
strains, i.e. 



n s (r s ) x u(r s ) x n s (r s ) = 0, 



(6) 



that can be easily deduced from the definition of the tan- 
gent component of the strain tensor. In (@), (|^) ns(rs) 
is the surface normal at a given point. 

Thus, the problem of equilibrium strain distribution 
in the definite-shaped sample below the phase transition 
temperature should be reformulated as follows: equations 
(|l|),(||) should be satisfied inside the sample volume for 
r 6 V\S, and equations (j^), are fulfilled at the sample 
surface, r = rg € S. 

It is obvious that for the symmetry-governed tran- 
sition, onset of homogeneous strain inevitably breaks 
the condition (||). The only exceptions being the cases 
of a rather thin disk- or plate-shaped sample with the 
platep-surface coinciding with the homogeneous shear 
plandiM 



III. "STRAY" ENERGY 

The problem of the equilibrium DS in AFM can be for- 
mally solved in analogy with the well known problem of 
the domain distribution in ferromagnets; namely, we as- 
sume that equilibrium (and in general, inhomogeneous) 
shear strain of the sample that complies with the condi- 
tion m) consists of two parts: 



Wshear(r) = u ms (r) + w ch (r), 



(7) 



where the first term, u ms (r), corresponds to the 
eigen spontaneous magnetostriction-induced strain de- 
fined from the free energy minimum (Q), and the second 



term, u c h(r) is an additional strain field produced by the 
"incompatibility clastic charges" localized at the crystal 
surface with the density: 

e c ias(r) = -n s x u ms (r s ) x n s 5/[n s (r - r s )], (8) 

where S is the Dirac's delta-function, prime means deriva- 
tive vs argument. Function tt c h(r) can be found explicitly 
from equation 



inc u ch (r) = e c i as (r) 
as follows (for the details see Ref. [n]) 

„ , s 1 f , e e ias(ri) - iTre e ias(ri) 
«ch(r) = — / dri 



(9) 



in 



v 



T — Tl 



1 f JC (n,r-rg) ~ 

— / dS- U r 5 

4tt J s \t-ts \ 6 

1 



in 



dn,U(r s ), 



(10) 



where 



U s (r) = (r) + n s (g) n s TrM ms (r) 
-n s (8 (u ms (r)n s ) - (u ms (r)n s ) ® n s . 

and dQ r is an increment of solid angle at which the sur- 
face point rs is seen from the given point r. 
Strain w c h(r) produces the "twinning" stress 



o-(r) = cu ch (r), 



(11) 



completely analogous to demagnetization field in ferro- 
magnets. In ( pd[ ) c is the tensor of the elastic mod- 
ula. Equilibrium strain distribution conditioned by non- 
deformed surface can thus be found from equation: 



du(r) 



cuch(r), 



(12) 



that combines ( |iT| ) and (Q) . It is easy to see that ( [l2| ) can 
be also deduced from variation of free energy functional 
with respect to the strain tensor components that satisfy 
compatibility conditions ([j]) inside the sampleED: 



$= / drFei(r) -F stray , 
Jv 



(13) 



where 



1 



-Fstray = ~ I dv / dS 



-in 



(n,r-r g ) 
l(r-rs) I* 



u(r)5f7(r s ) (14) 



is a twinning (stray) energy that was not considered in 
Ref. |. 

The above- developed approach makes it possible to 
solve the complicated problem of the domain distribu- 
tion in the nonlinear media under the action of external 
field in terms of strain tensor which in contrast to the 
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shift vector is observable and symmetrically related with 
an order parameter. 

To continue the analogy with the problem of ferromag- 
netic domains, we should emphasize some properties of 
the twinning (stray) fields in thermoelastics. First, as it 
is clearly seen from (|l0|), "incompatibility charges" can 
produce homogeneous additional strain u c h and thus, a 
macrostress in the case of appropriate form of the sam- 
ple. What is important, the energy contribution resulting 
from the macrostress is proportional to the sample vol- 
ume V and thus is principal, regardless of sample size. 
On the contrary to ferromagnets, where demagnetiza- 
tion charges are the scalars, in the case of elasticity, the 
charges (see (||)) possess the tensor characteristics, so, 
the only obvious case for homogeneous twinning stress is 
the thin plate. This case will be considered in details be- 
low. The other interesting feature of the additional strain 
field ( |l0| ) is that it is scaling-invariant. In other words, 
additional strain distribution inside the sample depends 
only upon the angle at which the surface is seen from the 
given point. So, isomorphic transformation of the sample 
does not change the additional strain distribution. 

From the above considerations we can make rather 
general conclusion: in the case of temperature-induced 
phase transition the macroscopic symmetry of low- and 
high temperature phases is the same. If the transition 
is symmetry-breaking (at the microscopic scale), i.e. the 
microscopic order parameter is conjugated with noniso- 
morphous striction and produces the strains that locally 
reduce the symmetry of the crystal lattice, then, the 
macroscopic symmetry is restored due to the onset of 
the DS with the differently oriented strain tensors. 



IV. EQUILIBIRUM DOMAIN STRUCTURE 



TP — T 2 rl 

r stray — l^/ Li 



Indeed, expressions (|13|) and (|14|) show that the shape 
of the sample essentially affects the strain distribution 
below the phase transition point. The simplest way to 
elucidate this point is to consider the case of thin plate 
for which the additional strains and twinning stresses are 
homogeneous. Suppose, the plate surface is rig-oriented 
with respect to crystal axes in high-symmetry phase, the 
transverse dimensions L of the plate are much greater 
than its thickness d. The main contribution to the stray 
energy (O) arises from the averaged magnetostriction: 



= t7 / dV ti™.s 
v Jv 



(r). 



(15) 



Substituting (|T|) into © and ( |TT| ) , one obtains a simple 
expression for the "charge" -induced macrostresses: 



CTmacro = c(U(r S )). 



(16) 



Altimately, from (|Rj) and (^4() stray energy is obtained 
as follows: 



|(u ms )c(?i ms ) + n s <g> n s Tr(u ms ) 
n s ® ((u mB }n s ) - ((u ms )n s ) <8> n s } . (17) 



If we restrict ourselves with the shear strains u ms = 
Mshcar, as it was assumed above (see (0) and (R)), then 



the stray energy takes a form 



Astray = L 2 d(u ms )c(u ms ) > 0, 



(18) 



that explicitly reveals it positivity for any choice of shear 
strain. In analogy with the case of ferromagnets, in the 
absence of external field the stray energy -F s tray can only 
be reduced by zeroing of the average strains. Expression 
( |l8| ) coincides with the expression (2) of Ref. || written 
previously from the phenomenologic considerations only. 

The value of F str ay (^8|) obviously depends upon the 
shape of the sample, in our case, from the orientation of 
plate surface with respect to crystal axes. If the plate 
normal rig is directed along the principal symmetry axis 
C n of the sample (n=3, 4, or 6) then, the macroscopic 
symmetry of the sample must be restored below the phase 
transition point and all the types of the domains should 
be equally represented. 

The above calculations convincingly show that appear- 
ance of the homogenous deviatoric (shear) strains in the 
finite-size sample can give rise to a considerable energy 
increase and thus are non- advantageous. There is a close 
analogy between the appearance of the long-range elastic 
fields in thermoelastics and dipole fields in ferromagnets 
(a similarity between the equations for electro-magnetic 
fields in substances and equations of elasticity theory was 
long ago noticed by R. de Wit in Ref. [r]) . The case of the 
weak ferromagnets needs the special treatment because 
the DS in corresponding compounds can be formed due 
to competition of demagnetization and twinning factors. 
This problem is out of scope of this paper. 

The general expression (|l^) can also be used for the 
description of the martensitic phase transitions with 
the only corrections concerning adjustable habitus (zero- 
shift) plane instead of fixed sample surface. 

Another important question that arises while analyz- 
ing the expression (|l7) is whether onset of twin (do- 
main) structure is indeed thcrmodynamically advanta- 
geous. The matter is that the inhomogeneous strain dis- 
tribution gives rise to an increase of the full free energy 
due to the contribution of the domain walls. In order to 
show that the energy gain ( [l8| ) because of the twinning is 
greater than the loss from the domain walls contribution, 
let us follow the method applied for the ferromagnets and 
consider the DS consisting of two alternating types of the 
domains, characterized with the strain tensors u\ and 112, 
the period of the structure g?ds is much less then the plate 
thickness, dns "C d. We suppose that the domains are 
coherently conjugated with each other and go out to the 
surface with the same strains that in the bulk. 

Standard calculations based on the formulae ( |l3| ) and 
( |l4| ) give rise to the following pantribution into free en- 
ergy from the "elastic charges"!^ 
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F, 



ch 



i 2 rfDsCl6("l - U 2 )c(ui - U 2 ) COS 



(19) 



where $ is the angle between the plate surface and do- 
main (twin) interface, £i(= 1 — £2) is the volume fraction 
of the domain of the first type. 

Domain walls (interfaces) also contribute into the free 
energy, corresponding expression is: 



Few — L °dw 



sfos ' 



(20) 



where odw is the surface energy of a single domain wall. 

Comparison of equations ( p7| ) , ( |l8| ) and ( |l9| ) shows that 
the energy of the domain walls compete with the contri- 
bution from the "elastic charges" and thus defines an 
optimal DS period: 



,opt _ 

"ds — 



\J 66 (ui - u 2 )£(mi - 112) cosd 



F e \ 



(21) 



So, contribution from the inhomogeneous part of 
strains and short-range periodic "charge" distribution is 
proportional to the ratio c£ds /d <C 1 and is much less than 
the stray energy (17) that is proportional to the sample 
volume. This means that onset of the DS in thermoelas- 
tics crucially diminishes the free energy from the "elas- 
tic charges" localized on the sample surface, the energy 
increase related with the inhomogeneous strain distribu- 
tion is much smaller and cannot compensate the energy 
related with homogezuzation of the domain structure for 
rather large sampleso. 



V. CONCRETE EXAMPLES 

To elucidate the consequences of the stray energy, let 
us consider a case of plate-shaped easy-plane pure AFM 
with the plane normal directed along the principal crystal 
axis. A good example of such an AFM is given by C0CI2 
(symmetry group is D^d) and low-doped YE^CuaOg+K 
at x < 0.3 (symmetry group is D^h) single crystals. Cor- 
responding contribution into stray energy ([17]) is: 

Astray = £ 2 d j^Cii [(u xx ) 2 + (u yv ) 2 ] 

+Cl2(u xx )(u yy ) + 2c 66 ( 1l X y ) 2 }, (22) 

where z-axis is supposed to be directed along the plate 
normal. Full free energy with the account of the magnetic 
subsystem in the external magnetic field H (neglecting 
the small demagnetization effects) could be written as 
follows 



F = L d 



k 



6[#a 



Mo 
2H K 



(m k f 



+UfcA mc lfc ® Ifc + -UkCUk] f + Astray, 



(23) 



where E an is the magnetic anisotropy energy, Mq is 
saturation magnetization, 2He is the value of spin-flip 
field of exchange nature, 1& and Uk are the vector of 
AFM and spontaneous magnetostriction, respectively, in 
the fc-th domain (k—1, 2, 3 for C0CI2 and k=l, 2 for 
YBa2Cu30e+a;), 6 is the volume fraction of the A:-th 

domain, A me is the 4-rank tensors of magnetoelastic co- 
efficient, the strain tensor is averaged over the sample 
volume as follows: < u >= £kp>k- 

Minimization of the energy ([23|) with respect to com- 
ponents of AFM vectors, strain tensor and also to domain 
fraction 6 provesQ that due to redistribution of the do- 
mains, the effective magnetic field inside sample is zero 
until the external field attains the value 



H \=Hi 



MD 



AcffMo 



HeMq 

Ceff 



(24) 



where A c ff and c e g are certain combinations of magne- 
toelastic and elastic constants that should be calculated 
with the account of concrete symmetry of the crystal. 
The value Hmd can be associated with the field of mon- 
odomenization of the sampleO. 

Below the filed of monodomenization the averaged 
shear strain 



\Wshca17 



AefFM 2 f H 



Ceff 



H 



MD 



= a 



H 



H 



MD 



(25) 



follows the quadratic field dependence normalized to the 
factor of monodomenization field value. In (25) uq de- 
scribes the absolute value of the local spontaneous shear. 
Thus, the low-field dependence of any macroscopic prop- 
erty which depends upon the relative fraction of the elas- 
tic domains should follow the law of the "corresponding 
states"! 2 -!!, i.e., should coincide for different substances in 
the reduced coordinates H/Hu,t> f° r the case H < Hmb- 



VI. CONCLUSIONS 

The main results of the paper can be formulated as 
follows: 

1. The equilibrium DS that arises in the course of 
thermoelastic phase transition results from the 
finite-size effects closely related with the proper- 
ties of sample surface. The absolutely rigid sur- 
face produces the imaginary "incompatibility elas- 
tic charges" which long-range field contributes into 
the energy of the sample, corresponding contribu- 
tion being proportional to the sample volume and 
thus provides the twinning effect. 

2. Contribution from inhomogeneous distribution of 
the strain below the phase transition point is pro- 
portional to the DS period and in the case of large 
samples is vanishingly small. 
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3. Elastic stray energy is the reason of the onset of 
equilibrium DS in pure AFM. The domain distri- 
bution can be reversibly regulated by the external 
magnetic field. 

4. Low-field dependence of macroscopic parameters 
should be the same for different samples if com- 
pared in the reduced coordinates H/Hmd (the law 
of the "corresponding states" ) . 
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